Taxonomic distribution of metabolic functions in bacteria associated with Trichodesmium consortia.

IMPORTANCE: Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions.

Oligotrophic waters of the Northwest Atlantic support taxonomically diverse diatom communities that are distinct from coastal waters.

Diatoms are important components of the marine food web and one of the most species-rich groups of phytoplankton. The diversity and composition of diatoms in eutrophic nearshore habitats have been well documented due to the outsized influence of diatoms on coastal ecosystem functioning. In contrast, patterns of both diatom diversity and community composition in offshore oligotrophic regions where diatom biomass is low have been poorly resolved. To compare the diatom diversity and community composition in oligotrophic and eutrophic waters, diatom communities were sampled along a 1,250 km transect from the oligotrophic Sargasso Sea to the coastal waters of the northeast US shelf. Diatom community composition was determined by amplifying and sequencing the 18S rDNA V4 region. Of the 301 amplicon sequence variants (ASVs) identified along the transect, the majority (70%) were sampled exclusively from oligotrophic waters of the Gulf Stream and Sargasso Sea and included the genera Bacteriastrum, Haslea, Hemiaulus, Pseudo-nitzschia, and Nitzschia. Diatom ASV richness did not vary along the transect, indicating that the oligotrophic Sargasso Sea and Gulf Stream are occupied by a diverse diatom community. Although ASV richness was similar between oligotrophic and coastal waters, diatom community composition in these regions differed significantly and was correlated with temperature and phosphate, two environmental variables known to influence diatom metabolism and geographic distribution. In sum, oligotrophic waters of the western North Atlantic harbor diverse diatom assemblages that are distinct from coastal regions, and these open ocean diatoms warrant additional study, as they may play critical roles in oligotrophic ecosystems.

Microbial metabolites in the marine carbon cycle.

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth.

Complex marine microbial communities partition metabolism of scarce resources over the diel cycle.

Complex assemblages of microbes in the surface ocean are responsible for approximately half of global carbon fixation. The persistence of high taxonomic diversity despite competition for a small suite of relatively homogeneously distributed nutrients, that is, 'the paradox of the plankton', represents a long-standing challenge for ecological theory. Here we find evidence consistent with temporal niche partitioning of nitrogen assimilation processes over a diel cycle in the North Pacific Subtropical Gyre. We jointly analysed transcript abundances, lipids and metabolites and discovered that a small number of diel archetypes can explain pervasive periodic dynamics. Metabolic pathway analysis of identified diel signals revealed asynchronous timing in the transcription of nitrogen uptake and assimilation genes among different microbial groups-cyanobacteria, heterotrophic bacteria and eukaryotes. This temporal niche partitioning of nitrogen uptake emerged despite synchronous transcription of photosynthesis and central carbon metabolism genes and associated macromolecular abundances. Temporal niche partitioning may be a mechanism by which microorganisms in the open ocean mitigate competition for scarce resources, supporting community coexistence.

Aureococcus anophagefferens (Pelagophyceae) genomes improve evaluation of nutrient acquisition strategies involved in brown tide dynamics.

The pelagophyte Aureococcus anophagefferens causes harmful brown tide blooms in marine embayments on three continents. Aureococcus anophagefferens was the first harmful algal bloom species to have its genome sequenced, an advance that evidenced genes important for adaptation to environmental conditions that prevail during brown tides. To expand the genomic tools available for this species, genomes for four strains were assembled, including three newly sequenced strains and one assembled from publicly available data. These genomes ranged from 57.11 to 73.62 Mb, encoding 13,191-17,404 potential proteins. All strains shared ~90% of their encoded proteins as determined by homology searches and shared most functional orthologs as determined by KEGG, although each strain also possessed coding sequences with unique functions. Like the original reference genome, the genomes assembled in this study possessed genes hypothesized to be important in bloom proliferation, including genes involved in organic compound metabolism and growth at low light. Cross-strain informatics and culture experiments suggest that the utilization of purines is a potentially important source of organic nitrogen for brown tides. Analyses of metatranscriptomes from a brown tide event demonstrated that use of a single genome yielded a lower read mapping percentage (~30% of library reads) as compared to a database generated from all available genomes (~43%), suggesting novel information about bloom ecology can be gained from expanding genomic space. This work demonstrates the continued need to sequence ecologically relevant algae to understand the genomic potential and their ecology in the environment.

Microbial community transcriptional patterns vary in response to mesoscale forcing in the North Pacific Subtropical Gyre.

The physical and biological dynamics that influence phytoplankton communities in the oligotrophic ocean are complex, changing across broad temporal and spatial scales. Eukaryotic phytoplankton (e.g., diatoms), despite their relatively low abundance in oligotrophic waters, are responsible for a large component of the organic matter flux to the ocean interior. Mesoscale eddies can impact both microbial community structure and function, enhancing primary production and carbon export, but the mechanisms that underpin these dynamics are still poorly understood. Here, mesoscale eddy influences on the taxonomic diversity and expressed functional profiles of surface communities of microeukaryotes and particle-associated heterotrophic bacteria from the North Pacific Subtropical Gyre were assessed over 2 years (spring 2016 and summer 2017). The taxonomic diversity of the microeukaryotes significantly differed by eddy polarity (cyclonic versus anticyclonic) and between sampling seasons/years and was significantly correlated with the taxonomic diversity of particle-associated heterotrophic bacteria. The expressed functional profile of these taxonomically distinct microeukaryotes varied consistently as a function of eddy polarity, with cyclones having a different expression pattern than anticyclones, and between sampling seasons/years. These data suggest that mesoscale forcing, and associated changes in biogeochemistry, could drive specific physiological responses in the resident microeukaryote community, independent of species composition.

Metatranscriptome Library Preparation Influences Analyses of Viral Community Activity During a Brown Tide Bloom.

There is growing interest in the use of metatranscriptomics to study virus community dynamics. We used RNA samples collected from harmful brown tides caused by the eukaryotic alga Aureococcus anophagefferens within New York (United States) estuaries and in the process observed how preprocessing of libraries by either selection for polyadenylation or reduction in ribosomal RNA (rRNA) influenced virus community analyses. As expected, more reads mapped to the A. anophagefferens genome in polyadenylation-selected libraries compared to the rRNA-reduced libraries, with reads mapped in each sample correlating to one another regardless of preprocessing of libraries. Yet, this trend was not seen for reads mapping to the Aureococcus anophagefferens Virus (AaV), where significantly more reads (approximately two orders of magnitude) were mapped to the AaV genome in the rRNA-reduced libraries. In the rRNA-reduced libraries, there was a strong and significant correlation between reads mappings to AaV and A. anophagefferens. Overall, polyadenylation-selected libraries produced fewer viral contigs, fewer reads mapped to viral contigs, and different proportions across viral realms and families, compared to their rRNA-reduced pairs. This study provides evidence that libraries generated by rRNA reduction and not selected for polyadenylation are more appropriate for quantitative characterization of viral communities in aquatic ecosystems by metatranscriptomics.

Combined pigment and metatranscriptomic analysis reveals highly synchronized diel patterns of phenotypic light response across domains in the open oligotrophic ocean.

Sunlight is the most important environmental control on diel fluctuations in phytoplankton activity, and understanding diel microbial processes is essential to the study of oceanic biogeochemical cycles. Yet, little is known about the in situ temporal dynamics of phytoplankton metabolic activities and their coordination across different populations. We investigated diel orchestration of phytoplankton activity in photosynthesis, photoacclimation, and photoprotection by analyzing pigment and quinone distributions in combination with metatranscriptomes in surface waters of the North Pacific Subtropical Gyre (NPSG). We found diel cycles in pigment abundances resulting from the balance of their synthesis and consumption. These dynamics suggest that night represents a metabolic recovery phase, refilling cellular pigment stores, while photosystems are remodeled towards photoprotection during daytime. Transcript levels of genes involved in photosynthesis and pigment metabolism had synchronized diel expression patterns among all taxa, reflecting the driving force light imparts upon photosynthetic organisms in the ocean, while other environmental factors drive niche differentiation. For instance, observed decoupling of diel oscillations in transcripts and related pigments indicates that pigment abundances are modulated by environmental factors extending beyond gene expression/regulation reinforcing the need to combine metatranscriptomics with proteomics and metabolomics to fully understand the timing of these critical processes in situ.

Transcriptional response of Emiliania huxleyi under changing nutrient environments in the North Pacific Subtropical Gyre.

The widespread coccolithophore Emiliania huxleyi is an abundant oceanic phytoplankton, impacting the global cycling of carbon through both photosynthesis and calcification. Here, we examined the transcriptional responses of populations of E. huxleyi in the North Pacific Subtropical Gyre to shifts in the nutrient environment. Using a metatranscriptomic approach, nutrient-amended microcosm studies were used to track the global metabolism of E. huxleyi. The addition of nitrate led to significant changes in transcript abundance for gene pathways involved in nitrogen and phosphorus metabolism, with a decrease in the abundance of genes involved in the acquisition of nitrogen (e.g. N-transporters) and an increase in the abundance of genes associated with phosphate acquisition (e.g. phosphatases). Simultaneously, after the addition of nitrate, genes associated with calcification and genes unique to the diploid life stages of E. huxleyi significantly increased. These results suggest that nitrogen is a major driver of the physiological ecology of E. huxleyi in this system and further suggest that the addition of nitrate drives shifts in the dominant life-stage of the population. Together, these results underscore the importance of phenotypic plasticity to the success of E. huxleyi, a characteristic that likely underpins its ability to thrive across a variety of marine environments.

Progress and promise of omics for predicting the impacts of climate change on harmful algal blooms.

Climate change is predicted to increase the severity and prevalence of harmful algal blooms (HABs). In the past twenty years, omics techniques such as genomics, transcriptomics, proteomics and metabolomics have transformed that data landscape of many fields including the study of HABs. Advances in technology have facilitated the creation of many publicly available omics datasets that are complementary and shed new light on the mechanisms of HAB formation and toxin production. Genomics have been used to reveal differences in toxicity and nutritional requirements, while transcriptomics and proteomics have been used to explore HAB species responses to environmental stressors, and metabolomics can reveal mechanisms of allelopathy and toxicity. In this review, we explore how omics data may be leveraged to improve predictions of how climate change will impact HAB dynamics. We also highlight important gaps in our knowledge of HAB prediction, which include swimming behaviors, microbial interactions and evolution that can be addressed by future studies with omics tools. Lastly, we discuss approaches to incorporate current omics datasets into predictive numerical models that may enhance HAB prediction in a changing world. With the ever-increasing omics databases, leveraging these data for understanding climate-driven HAB dynamics will be increasingly powerful.

Transcriptional patterns of Emiliania huxleyi in the North Pacific Subtropical Gyre reveal the daily rhythms of its metabolic potential.

Emiliania huxleyi is a calcifying haptophyte, contributing to both the organic and inorganic marine carbon cycles. In marine ecosystems, light is a major driver of phytoplankton physiology and ultimately carbon flow through the ecosystem. Here, we analysed a Lagrangian time-series of metatranscriptomes collected in the North Pacific Subtropical Gyre (NPSG) to examine how in situ populations of E. huxleyi modulate gene expression over day-night transitions. Many E. huxleyi contigs had a diel expression pattern, with 61% of contigs clustering into modules with statistically significant diel periodicity. Contigs involved in processes that build up energy stores, like carbon fixation and lipid synthesis, peaked around dawn. In contrast, contigs involved in processes that released energy stores, like respiration and lipid degradation, peaked mid-day and towards dusk. These patterns suggest an orchestrated cycle of building, then consuming energy stores in E. huxleyi populations in the NPSG. Selected contigs related to the cell cycle also exhibited significant diel periodicity consistent with phased modulations of division observed in culture. Overall, these patterns of gene expression suggest a daily metabolic cascade that could contribute to both organic and inorganic carbon flow in this nutrient depleted ecosystem.

Kilauea lava fuels phytoplankton bloom in the North Pacific Ocean.

From June to August 2018, the eruption of Kilauea volcano on the island of Hawai'i injected millions of cubic meters of molten lava into the nutrient-poor waters of the North Pacific Subtropical Gyre. The lava-impacted seawater was characterized by high concentrations of metals and nutrients that stimulated phytoplankton growth, resulting in an extensive plume of chlorophyll a that was detectable by satellite. Chemical and molecular evidence revealed that this biological response hinged on unexpectedly high concentrations of nitrate, despite the negligible quantities of nitrogen in basaltic lava. We hypothesize that the high nitrate was caused by buoyant plumes of nutrient-rich deep waters created by the substantial input of lava into the ocean. This large-scale ocean fertilization was therefore a unique perturbation event that revealed how marine ecosystems respond to exogenous inputs of nutrients.

Transcriptional and Proteomic Choreography Under Phosphorus Deficiency and Re-supply in the N2 Fixing Cyanobacterium Trichodesmium erythraeum.

The N2 fixing cyanobacterium Trichodesmium is a critically important organism in oligotrophic marine ecosystems, supplying "new" nitrogen (N) to the otherwise N-poor tropical and subtropical regions where it occurs. Low concentrations of phosphorus (P) in these regions can constrain Trichodesmium distribution and N2 fixation rates. Physiological characterization of a single species in a mixed community can be challenging, and 'omic approaches are increasingly important tools for tracking nutritional physiology in a taxon-specific manner. As such, studies examining the dynamics of gene and protein markers of physiology (e.g., nutrient stress) are critical for the application and interpretation of such 'omic data in situ. Here we leveraged combined transcriptomics, proteomics, and enzyme activity assays to track the physiological response of Trichodesmium erythraeum IMS101 to P deficiency and subsequent P re-supply over 72 h of sampling. P deficiency resulted in differential gene expression, protein abundance, and enzyme activity that highlighted a synchronous shift in P physiology with increases in the transcripts and corresponding proteins for hydrolyzing organic phosphorus, taking up phosphate with higher affinity, and modulating intracellular P demand. After P deficiency was alleviated, gene expression of these biomarkers was reduced to replete levels within 4 h of P amendment. A number of these gene biomarkers were adjacent to putative pho boxes and their expression patterns were similar to a sphR response regulator. Protein products of the P deficiency biomarkers were slow to decline, with 84% of the original P deficient protein set still significantly differentially expressed after 72 h. Alkaline phosphatase activity tracked with proteins for this enzyme. With the rapid turnover time of transcripts, they appear to be good biomarkers of a P stress phenotype, whereas proteins, with a slower turnover time, may better reflect cellular activities. These results highlight the importance of validating and pairing transcriptome and proteome data that can be applied to physiological studies of key species in situ.

Transcriptional Shifts Highlight the Role of Nutrients in Harmful Brown Tide Dynamics.

Harmful algal blooms (HABs) threaten ecosystems and human health worldwide. Controlling nitrogen inputs to coastal waters is a common HAB management strategy, as nutrient concentrations often suggest coastal blooms are nitrogen-limited. However, defining best nutrient management practices is a long-standing challenge: in part, because of difficulties in directly tracking the nutritional physiology of harmful species in mixed communities. Using metatranscriptome sequencing and incubation experiments, we addressed this challenge by assaying the in situ physiological ecology of the ecosystem destructive alga, Aureococcus anophagefferens. Here we show that gene markers of phosphorus deficiency were expressed in situ, and modulated by the enrichment of phosphorus, which was consistent with the observed growth rate responses. These data demonstrate the importance of phosphorus in controlling brown-tide dynamics, suggesting that phosphorus, in addition to nitrogen, should be evaluated in the management and mitigation of these blooms. Given that nutrient concentrations alone were suggestive of a nitrogen-limited ecosystem, this study underscores the value of directly assaying harmful algae in situ for the development of management strategies.

Non-linear Physiology and Gene Expression Responses of Harmful Alga Heterosigma akashiwo to Rising CO2.

Heterosigma akashiwo is a raphidophyte known for forming ichthyotoxic blooms. In order to predict the potential impacts of rising CO2 on H. akashiwo it is necessary to understand the factors influencing growth rates over a range of CO2 concentrations. Here we examined the physiology and gene expression response of H. akashiwo to concentrations from 200 to 1000ppm CO2. Growth rate data were combined from this and previous studies and fit with a CO2 limitation-inhibition model that revealed an apparent growth optimum around 600-800ppm CO2. Physiological changes included a significant increase in C:N ratio at ∼800ppm CO2 and a significant decrease in hydrogen peroxide concentration at ∼1000ppm. Whole transcriptome sequencing of H. akashiwo revealed sharp distinctions in metabolic pathway gene expression between ∼600 and ∼800ppm CO2. Hierarchical clustering by co-expression identified groups of genes with significant correlations to CO2 and growth rate. Genes with significant differential expression with CO2 included carbon concentrating mechanism genes such as beta-carbonic anhydrases and a bicarbonate transporter, which may underpin shifts in physiology. Genes involved in cell motility were significantly changed by both elevated CO2 and growth rate, suggesting that future ocean conditions could modify swimming behavior in this species.

Periodic and coordinated gene expression between a diazotroph and its diatom host.

In the surface ocean, light fuels photosynthetic carbon fixation of phytoplankton, playing a critical role in ecosystem processes including carbon export to the deep sea. In oligotrophic oceans, diatom-diazotroph associations (DDAs) play a keystone role in ecosystem function because diazotrophs can provide otherwise scarce biologically available nitrogen to the diatom host, fueling growth and subsequent carbon sequestration. Despite their importance, relatively little is known about the nature of these associations in situ. Here we used metatranscriptomic sequencing of surface samples from the North Pacific Subtropical Gyre (NPSG) to reconstruct patterns of gene expression for the diazotrophic symbiont Richelia and we examined how these patterns were integrated with those of the diatom host over day-night transitions. Richelia exhibited significant diel signals for genes related to photosynthesis, N2 fixation, and resource acquisition, among other processes. N2 fixation genes were significantly co-expressed with host nitrogen uptake and metabolism, as well as potential genes involved in carbon transport, which may underpin the exchange of nitrogen and carbon within this association. Patterns of expression suggested cell division was integrated between the host and symbiont across the diel cycle. Collectively these data suggest that symbiont-host physiological ecology is strongly interconnected in the NPSG.

Shifting metabolic priorities among key protistan taxa within and below the euphotic zone.

A metatranscriptome study targeting the protistan community was conducted off the coast of Southern California, at the San Pedro Ocean Time-series station at the surface, 150 m (oxycline), and 890 m to link putative metabolic patterns to distinct protistan lineages. Comparison of relative transcript abundances revealed depth-related shifts in the nutritional modes of key taxonomic groups. Eukaryotic gene expression in the sunlit surface environment was dominated by phototrophs, such as diatoms and chlorophytes, and high abundances of transcripts associated with synthesis pathways (e.g., photosynthesis, carbon fixation, fatty acid synthesis). Sub-euphotic depths (150 and 890 m) exhibited strong contributions from dinoflagellates and ciliates, and were characterized by transcripts relating to digestion or intracellular nutrient recycling (e.g., breakdown of fatty acids and V-type ATPases). These transcriptional patterns underlie the distinct nutritional modes of ecologically important protistan lineages that drive marine food webs, and provide a framework to investigate trophic dynamics across diverse protistan communities.

Transcriptional patterns identify resource controls on the diazotroph Trichodesmium in the Atlantic and Pacific oceans.

The N2-fixing cyanobacterium Trichodesmium is intensely studied because of the control this organism exerts over the cycling of carbon and nitrogen in the low nutrient ocean gyres. Although iron (Fe) and phosphorus (P) bioavailability are thought to be major drivers of Trichodesmium distributions and activities, identifying resource controls on Trichodesmium is challenging, as Fe and P are often organically complexed and their bioavailability to a single species in a mixed community is difficult to constrain. Further, Fe and P geochemistries are linked through the activities of metalloenzymes, such as the alkaline phosphatases (APs) PhoX and PhoA, which are used by microbes to access dissolved organic P (DOP). Here we identified significant correlations between Trichodesmium-specific transcriptional patterns in the North Atlantic (NASG) and North Pacific Subtropical Gyres (NPSG) and patterns in Fe and P biogeochemistry, with the relative enrichment of Fe stress markers in the NPSG, and P stress markers in the NASG. We also observed the differential enrichment of Fe-requiring PhoX transcripts in the NASG and Fe-insensitive PhoA transcripts in the NPSG, suggesting that metalloenzyme switching may be used to mitigate Fe limitation of DOP metabolism in Trichodesmium. This trait may underpin Trichodesmium success across disparate ecosystems.

Coordinated gene expression between Trichodesmium and its microbiome over day-night cycles in the North Pacific Subtropical Gyre.

Trichodesmium is a widespread, N2 fixing marine cyanobacterium that drives inputs of newly fixed nitrogen and carbon into the oligotrophic ecosystems where it occurs. Colonies of Trichodesmium ubiquitously occur with heterotrophic bacteria that make up a diverse microbiome, and interactions within this Trichodesmium holobiont could influence the fate of fixed carbon and nitrogen. Metatranscriptome sequencing was performed on Trichodesmium colonies collected during high-frequency Lagrangian sampling in the North Pacific Subtropical Gyre (NPSG) to identify possible interactions between the Trichodesmium host and microbiome over day-night cycles. Here we show significantly coordinated patterns of gene expression between host and microbiome, many of which had significant day-night periodicity. The functions of the co-expressed genes suggested a suite of interactions within the holobiont linked to key resources including nitrogen, carbon, and iron. Evidence of microbiome reliance on Trichodesmium-derived vitamin B12 was also detected in co-expression patterns, highlighting a dependency that could shape holobiont community structure. Collectively, these patterns of expression suggest that biotic interactions could influence colony cycling of resources like nitrogen and vitamin B12, and decouple activities, like N2 fixation, from typical abiotic drivers of Trichodesmium physiological ecology.

The impact of elevated CO2 on Prochlorococcus and microbial interactions with 'helper' bacterium Alteromonas.

Prochlorococcus is a globally important marine cyanobacterium that lacks the gene catalase and relies on 'helper' bacteria such as Alteromonas to remove reactive oxygen species. Increasing atmospheric CO2 decreases the need for carbon concentrating mechanisms and photorespiration in phytoplankton, potentially altering their metabolism and microbial interactions even when carbon is not limiting growth. Here, Prochlorococcus (VOL4, MIT9312) was co-cultured with Alteromonas (strain EZ55) under ambient (400 p.p.m.) and elevated CO2 (800 p.p.m.). Under elevated CO2, Prochlorococcus had a significantly longer lag phase and greater apparent die-offs after transfers suggesting an increase in oxidative stress. Whole-transcriptome analysis of Prochlorococcus revealed decreased expression of the carbon fixation operon, including carboxysome subunits, corresponding with significantly fewer carboxysome structures observed by electron microscopy. Prochlorococcus co-culture responsive gene 1 had significantly increased expression in elevated CO2, potentially indicating a shift in the microbial interaction. Transcriptome analysis of Alteromonas in co-culture with Prochlorococcus revealed decreased expression of the catalase gene, known to be critical in relieving oxidative stress in Prochlorococcus by removing hydrogen peroxide. The decrease in catalase gene expression was corroborated by a significant ~6-fold decrease in removal rates of hydrogen peroxide from co-cultures. These data suggest Prochlorococcus may be more vulnerable to oxidative stress under elevated CO2 in part from a decrease in ecosystem services provided by heterotrophs like Alteromonas. This work highlights the importance of considering microbial interactions in the context of a changing ocean.The ISME Journal advance online publication, 31 October 2017; doi:10.1038/ismej.2017.189.